Do all states undergo sudden death of entanglement at finite temperature? 
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In this paper we consider the decay of quantum entanglement, quantified by the concurrence, of 
a pair of two-level systems each of which is interacting with a reservoir at finite temperature T. 
For a broad class of initially entangled states, we demonstrate that the system always becomes 
disentangled in a finite time i.e. "entanglement sudden death" (ESD) occurs. This class includes 
all states which previously had been found to have long-lived entanglement in zero temperature 
reservoirs. Our general result is illustrated by an example. 



In the past few years there has been considerable in- 
terest in the properties of entangled quantum systems. 
Spurred on by the emergence of compelling applications 
in quantum information processing, useful methods by 
which the entanglement of quantum systems can be es- 
tablished and characterized have emerged. Perhaps the 
most impactful to date has been the simple procedure 
derived by Wootters 1 for quantifying entanglement for 
an arbitrary mixed state of a pair of two-level systems. 
This has provided a very useful tool for measurement of 
experimental quantum states [2] and is today commonly 
used in assessing the capabilities of emerging quantum 
technologies. Building on Wootter's work, recently Yu 
and Eberly [3J investigated the time evolution of entan- 
glement (quantified using the concurrence) of a bipartite 
qubit system undergoing various modes of decoherence. 
Remarkably, they found that, even when there is no in- 
teraction, (either directly or through a correlated envi- 
ronment), there are certain states whose entanglement 
decays exponentially with time, while for other closely 
related states, the entanglement vanishes complctly in a 
finite time. This "entanglement sudden death" (ESD) is 
an intriguing and potentially very important discovery. 
Since the first theoretical demonstration of ESD, further 
investigations of different systems have been made by 
various groups [H EJ El S HI El HOI HU C21 E3] ■ Extending 
Yu and Eberly's model by considering correlated reser- 
voirs and interactions [3 [TU1 [12 E], it was shown 
that entanglement may be created by spontaneous emis- 
sion (something which has been known for some time |17j 
in a different context) . The ESD has also been predicted 
for more complicated systems using other entanglement 
measures [T5J and an attempt to give a geometric 
interpretation for the phenomena of ESD has also been 
made [18]. Very recently, experimental studies have also 
been carried out to demostrate ESD, using carefully en- 
gineered interactions between systems and environments: 
Sudden death has been observed both in photons [H] and 
in atomic ensembles [19] . 

From the technological point of view, states which ex- 
hibit exponential decay of entanglement, and therefore 
retain some vestige of this all-important correlation for 
long periods, are of great interest. Thus it is important to 
identify precisely in what circumstances ESD will occur. 
In this paper, we consider qubits in finite temperature 



reservoirs: instead of the energy of the qubits being lost 
via spontaneous decay to the environment, now addition- 
ally the reservoirs can cause excitation of the qubits. For 
a broad class of mixed quantum states, which includes 
all of the states studied by Yu and Eberly and others 
in connection with this problem, we demonstrate that 
all states undergo sudden death of entanglement at finite 
temperature. 

As in [3J, we study a system of two qubits initially 
entangled and interacting with uncorrelated reservoirs. 
However, unlike |3J, in which the system is studied at 
T=0, we include the effects of heat in our system (Fig.[lJ. 
Here the dynamics of the density matrix p describing the 
two qubits is given by: 

^ = ±[H,p}+L 1 [p]+L 2 [p], (1) 

where [H, p] is the unitary part of the evolution (which we 
shall ignore as it has no effect on our study decoherence) . 
The Liouvillian of the ith qubit is given by: 

W = &^{[*iM] + [**-&&*+]} 

riV 

+ T {[^ + ,pa*_] + [&\p,&l\}, (2) 
where T is the spontaneous decay rate of the qubits, a l + = 




FIG. 1: Disentanglement by spontaneous emission of 
a two-qubit system in a heat bath. The reservoir is 
modelled by different harmonic oscillator modes. Each qubit, 
here depicted by a two-level atom, interacts with its reservoir. 
The only interaction between the qubits that ever exists is the 
one that leads to their entanglement at time t =0. Following 
this, however, the only interaction that remains is that with 
the corresponding reservoir. This leads to decoherence, which 
causes the qubits to disentangle. Here we include the effect 
of heat by studying the system at T / 0. 
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(|1> (0|)i, and at = (|0) (1|) 4 , where the index i (=1,2) 
denotes the qubits. The first term on the right hand 
side of ([2]) corresponds to the depopulation of the atoms, 
while the second term describes the decay caused by the 
finite temperature; and n is the mean occupation number 
of the reservoir (assumed to be the same for both) . 

We assume that our system is initially an"X-state" de- 
scribed by the following density matrix: 



z(t) 



z X, 



Pit) 



( a(t) w(t) \ 

b(t) z(t) 

z*{t) c(t) 

\ w*(t) d{t) ) 



(3) 



Such states are general enough to include states such as 
the Werner states, the maximally entangled mixed states 
(MEMS) [21], the Bell States, and what we will refer to as 
the pye states, studied in [3] and which will be described 
later. 

Substituting ^ into ([!]), the Master equation of our 
system, we obtain the following first order coupled dif- 
ferential equations: 

a(t) = T [-2{n + l)a(t) + b(t)n + c(t)n] , 
b(t) =T[(n+ l)a(i) - (2n + l)b(t) + nd(t)\ , 
c(t) = T [(n + l)a(t) - (2n + l)c(t) + nd(t)} , 
d(t) = T [(n + l)6(f) + (n + l)c(t) - 2nd(t)\ , 
z(t)=T[-(2n + l)z(t)], 
w(t) = r[-(2n+ l)w{t)] . 

(4) 

These may be solved to yield the following expressions: 



a(t) 



b(t) 



1 



c(t) 



{n(n+l) 



(5) 



where X 



-r(2n+i)t 



a =a(0), etc. 



(2n- 

+ [2(a - d )n 2 + (o - rf + 1)"-]^ 
+ [(2a + 2d - 1)« 2 + (3a + rf - l)n + a ]X 2 } 

= (2^TIF { ^ + 1) 

-[2(a + 2c + rf - 1)« 2 + (ao + 4c + 3d - 2)n 
+ (c + d -l)]^ 

-[(2a + 2d - l)n 2 + (3a + d - l)n + a ]X 2 } 
1 

~ (2n+ 1) : 

+ [2(a + 2c + rf - 1)« 2 + (3a + 4c + d a - 2)n]X 
-[(2a + 2d - 1)« 2 + (3a + rf - l)n + ao]^ 2 } , 



Using Wootters' formula [T] , the concurrence for a state 
of the form given by ^ is: 

C = 2max jo, \z(t)\ - y/a(t)d(t), \w(t)\ - y/b(t)c(t)} . 

(6) 

This implies that the disentanglement time will be the 
largest positive solutions of the following equations: 

\z(t)\ - y/a(t)d(t) = 0, Hi) I - y/b(t)c(t) - 0. (7) 

Multiplying both equations in (7| by the positive quanti- 
ties \z{t)\ + yj a{t)d{t) and \w(t) + yfb{t)c{t), respectively, 
yields: 

\z(t)\ 2 - a(t)d(t) = 0, \w(t)\ 2 - b(t)c(t) = 0. (8) 

Substituting from equation (JsJ) gives two quartic equa- 
tions in X, which we will use in the proof of our main 
result. 

The quantity X is the time-dependent parameter that 
we use to monitor the evolution of entanglement in the 
system. Notice that at t=0, X=l, and that at t=oo, 
X=0. Physically meaningful values for X are, therefore, 
between and 1 . Asymptotic decay of entanglement im- 
plies a solution at X=0. However, if the entanglement of 
the system decays in a finite time (ESD), the solution of 
^ must lie in the range < X < 1. 

Both equations in (JsJ) are polynomial equations in X 
and continuous. At X=0, these equations take the fol- 
lowing value: 



-n 2 {n 2 + 2n + l) 
(2n+ l) 4 ' 



(9) 



d(t) = 



1 



w(t) 



-(n+ l)[2n(ao - do) + (o - d + 1)]X 

+ [(2a + 2d - l)n 2 + (3a + d - l)h + a ]X 2 } 

= w X, 



Notice that since n is a positive quantity, ^ is negative. 
On the other hand, if we evaluate ^ for X=l, which 
corresponds to t=0, we obtain \z\ 2 — ad and \w\ 2 — be for 
the first and second equations, respectively. At least one 
of these quantities has to be positive if we assume that 
our systems are initially entangled, which is the case here. 
The fact that the quartic equations are continuous and 
have a negative value at X=0 and a positive one at X=l 
implies that they have at least one root in our interval of 
interest < X < 1. See Fig. [2j 

For finite n, and therefore for finite T, the constant 
term ^ is always finite and nonzero for h > 0. Hence, 
there is no X=0 solution; i.e., no asymptotic decay. 
Equation ([8| has at least one solution in the range 
< X < 1. implying that the entanglement falls to 
zero in a finite time. 
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As an example, let us consider a special case when 
w(t)=0. In this case, the only quartic equation that has 
to be satisfied for C=0 is the first one in ^ . This equa- 
tion can be easily solved with appropriate reparametriza- 
tion. The four solutions for X are given by: 

X = (r + s)± v/(r + s) 2 + < 2 , 

X = (r-s)± v/(r-s) 2 + < 2 . (10) 

where 

_ 1 + (go - d )(2n+ 1) 

T ~ 4 {(2n + 1) [a (n + 1) + d n] - n(n + 1)} ' 

= (2n + 1) V (1 - op - rfo) 2 + 4(z 2 - a do) 
S ~~ 4{(2n + 1) [a (n + 1) + d a n] - n(n + l)}' 




n(n + 1) 

{(2n + 1) [a (n+ 1) + d «-] = n(n+ 1)}' 



Further, following ([3|, let us consider states of the 
form: 



In this paper, we presented a proof that in two-qubit 
systems interacting with uncorrelated reservoirs and de- 
scribed by X-states, ESD always occurs at any finite tem- 
perature. Although X-states are quite general states, 
they are not the most general ones. Thus the next ques- 
tion to ask is: do all states exhibit ESD? In other words, 
for a state described by a density matrix, without any 
zero elements, is entanglement still lost in a finite time? 
This question is not straight forward. One has to find 
the equation for concurrence in that case, determine its 
order, study the properties of its coefficient, and from 
there maybe be able to comment on the nature of the 
roots, and, therefore, be able to predict what will hap- 
pen in the actual physical system. However, our results 
demonstrate that in any finite temperature reservoir, all 
states that have been shown to be long-lived in a zero- 
temperature bath, will undergo sudden death. 
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(11) 



In [3], Yu and Eberly have shown that for < a < 
|, the entanglement of this state is long-lived at zero 
temperature. Here, we will illustrate that as soon as h 
becomes finite, the range vanishes and there is no long- 
lived entanglement for any value of a. The plots in Fig. [3] 
illustrate this. As n increases, the range of values for a 
in which entanglement is long-lived diminishes until all 
values disappear. 




FIG. 2: Plot of F(X) vs X. This is the plot of the first 
quartic equation in Q for n =0.8, a =0.1, d = 0.05, z =0.3. 
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